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ORIGINAL ARTICLE 

Population Pharmacokinetics and Antimalarial 
Pharmacodynamics of Piperaquine in Patients With 
Plasmodium vivax Malaria in Thailand 

J Taming^ ^ P Thana\ AP Phyo^ KM Lwin^ W Hanpithakpong\ EA Ashley'■'^ NPJ Day^ F Nosten^ ^'^ and NJ White^ ^' 

Dihydroartemisinin-piperaquine is an effective drug in the treatment of Plasmodium falciparum and P. vivax malaria. 
The objective of this study was to evaluate the population pharmacokinetics and pharmacodynamics of piperaquine in patients 
with P. vivax malaria in Thailand after a standard regimen of dihydroartemisinin-piperaquine to determine whether residual 
piperaquine prevents or delays the emergence of P. vivax relapse. Sparse blood samples were collected from 116 patients. 
Piperaquine pharmacokinetics were described well by a three-compartment distribution model. Relapsing P. v/Vax malaria was 
accommodated by a constant baseline hazard (8.94 relapses/year) with the addition of a surge function in a fixed 3-week interval 
and a protective piperaquine effect. The results suggest that a large proportion of the first relapses were suppressed completely 
by residual piperaquine concentrations and that recurrences resulted mainly from emergence of the second or third relapse 
or from reinfection. This suggests a significant reduction in P. vivax morbidity when using dihydroartemisinin-piperaquine 
compared with other antimalarial drugs with shorter terminal postprophylactic effects. 
CPT Pharmacometrics Syst. Pharmacol. (2014) 3, e132; doi:10.1038/psp.2014.29; published online 27 August 2014 



Plasmodium vivax Is the most difficult of the five species 
of human malaria parasites to eliminate. It Is the most geo- 
graphically widespread human malaria parasite, causing an 
estimated 1 30-390 million cases every year.^ P. vivax Is trans- 
mitted within 95 countries In America, Africa, and Asla.^'^ In 
total, 2.85 billion people were estimated to be at risk of P. vivax 
transmission In 2009, especially In Central and Southeast 
Asia, where more than 90% of the world's P. wVaxtransmisslon 
occurs.^ Acute P. vivax malaria Is rarely life threatening, but it 
causes substantial morbidity In higher transmission areas as 
a result of multiple relapses from latent liver-stage parasites."* 
These repeated infections may result In life-threatenlngly 
severe anemia due to destruction of both infected and non- 
infected red blood cells. P. vivax Is an important cause of low 
birth weight.""^ 

Chloroqulne, a blood schlzontlcide, Is used In first-line 
treatment and prophylaxis of P. wVax malaria In most endemic 
areas. However, P. vivax resistance to chloroqulne, which was 
first observed In Papua New Guinea In 1989,'' Is increasing 
in many parts of the world." A combination of two drugs with 
different mechanisms of action Is preferred for the treatment 
of malaria to reduce the risk of the development of resistant 
parasite strains. The fixed combination of dihydroartemisinin- 
piperaquine (DHA-PQ) Is recommended by the World Health 
Organization as the first-line treatment of uncomplicated 
P. falciparum malaria*" and has replaced chloroqulne In the 
first-line treatment of resistant P. vivax malaria in Indonesia. 
DHA-PQ is currently available as a fixed-dose combination, 
administered once dally over 3 days. A meta-analysis pooled 
from 12 different studies in 6 countries between 2003 and 
2006 showed high efficacy against P. falciparum malaria 
(polymerase chain reaction-corrected cure rates of 98.7% at 



day 28) and good tolerance (4.8% total incidence of early 
vomiting) of DHA-PQ." DHA-PQ also has excellent efficacy 
against P. wVax malaria.*^"" 

Piperaquine Is distributed extensively Into tissues, result- 
ing in multiphasic distribution and a slow elimination from the 
systemic circulation (elimination half-life (f^,^) °f ~1 8-28 days). 
Similar pharmacokinetic properties have been reported in 
children but these pediatric studies suggest that small children 
need a higher weight-based dose to achieve exposure com- 
parable to that in older children.*'' Most piperaquine pharma- 
cokinetic studies have sampled venous plasma. However, 
capillary sampling is a desirable alternative for field pharma- 
cokinetic studies, especially in children. The population phar- 
macokinetics and pharmacodynamics of piperaquine have 
not been reported previously In patients with P. vivax malaria. 

In South East Asia and Oceania, P. vivax relapses at 
3-week intervals. Treatment with efficacious slowly eliminated 
antimalarials results in substantially delayed emergence of 
the relapse, but it has been unclear whether the first relapse 
is prevented (and so the second relapse Is observed) or 
delayed. The aim of this study was to characterize the popu- 
lation pharmacokinetics and antimalarial pharmacodynamics 
of piperaquine In the treatment of P. wVax malaria and to use 
this Information to determine whether piperaquine prevents 
or delays the first P. vivax relapse. 



RESULTS 

Safety and efficacy 

This study was conducted at Shoklo Malaria Research Unit 
clinics, Mae Sot, Thailand, an area of low and seasonal 
malaria transmission located at the Thalland-Myanmar 
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Table 1 Clinical outcomes and patient demographics of patients with Plasmodium vivax malaria enrolled in this study 





Efficacy study 


PK cohort' 


PK-PD cohort" 


Total number of patients 


241 


116 


62 


Total number of venous/capillary samples 


— 


435/356 


386/312 


Total dose of piperaquine base (mg/kg) 


30.8 (23.1-46.2) 


30.1 (27.2-36.1) 


30.1 (27.2-36.1) 


Continuous and categorical admission covariates 








Age (years) 


18 (1-57) 


23.5 (3-57) 


29 (11-56) 


Height (cm) 


154 (65-179) 


160 (88-176) 


160 (115-176) 


Body weight (kg) 


45 (7-74) 


50 (9-67) 


50 (23-67) 


Body temperature (°C) 


37 (36-41) 


37 (36-40.4) 


37 (36^0.1) 


Days of fever before admission 


2 (0-8) 


2 (1-8) 


2 (1-7) 


Hematocrit (%) 


37 (22-50) 


39 (22-50) 


39 (27-^9) 


Initial parasitemia (pi ') 


6,106 (80-70,255) 


7,185 (80-58,714) 


6,802 (80^5,166) 


Gametocytomia (yes/no) 


165/76 


65/51 


34/28 


Sex (male/female) 


169/72 


90/26 


25/37 


Treatment outcomes" 








Number of patients with no recurrent malaria 


120 


31 


31 


Number of patients with recurrent P. vivax malaria 


109 


75 


31 


Time to recurrent P. vivax malaria (days) 


56(14-65) 


56 (21-65) 


56 (35-65) 


Parasitemia at recurrent P. v/Vax malaria (pi ') 


352 (16-39,144) 


735 (16-39,144) 


735 (16-31,617) 


Number of patients with recurrent P. falciparum malaria 


10 


8 




Time to recurrent P. falciparum malaria (days) 


37 (28^9) 


40 (35-49) 




Number of patients with recurrent mixed malaria 


2 


2 




Time to recurrent mixed malaria (days) 


42 (35-49) 


42 (35-49) 





All values are given as median (range) unless otherwise specified. 
PK-PD, pharmacokinetics/pharmacodynamics. 

"Based on 62 patients enrolled for pharmacokinetic modeling and an additional 54 patients who provided a single plasma sample only at the time of recurrent malaria. 
"Based on 62 patients for whom both pharmacokinetic and pharmacodynamic data were available. =As confirmed by polymerase chain reaction genotyping. 



border. This pharmacokinetic study was nested into a larger 
efficacy trial in 500 patients with P. vivax malaria infections. 
Full clinical details have been published elsewhere and DHA- 
PQ was found to be an effective alternative treatment with 
a lower cumulative risk of recurrent P. vivax malaria at 9 
weeks compared with chloroquine (54.9 vs. 79.1%).^'' Two- 
hundred and fifty Karen and Burmese patients with P. vivax 
malaria infections were enrolled in the DHA-PQ-treated arm 
of the efficacy study (Table 1). All patients received DHA- 
PQ (Duo-cotecxin, Beijing Holley-Cotec Pharmaceuticals, 
China) tablets containing 40 mg of dihydroartemisinin and 
320 mg of piperaquine phosphate. A standard weight-based 
regimen (2.2mg/kg/day of dihydroartemisinin and 17.8mg/ 
kg/day of piperaquine phosphate) was administered with milk 
once daily for 3 days. The fixed combination of DHA-PQ was 
well tolerated with no serious adverse events reported dur- 
ing follow-up. In this population modeling, nine patients were 
excluded: two patients did not meet the inclusion criteria, five 
patients had an incomplete course of medication, one patient 
self-medicated with other drugs, and one patient withdrew 
consent. Overall, 45% (109/241) of patients presented with 
recurrent P. vivax malaria during the follow-up period at a 
median (range) of 56 (14-65) days (Table 1). 

Pharmacokinetics of piperaquine 

A total of 791 plasma concentrations (435 venous and 356 
capillary samples) from 116 patients were quantified and 
used in the pharmacokinetic modeling. Venous and capil- 
lary plasma concentrations were transformed into their 
natural logarithms and modeled simultaneously using 
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Figure 1 Final piperaquine population pharmacokinetic- 
pharmacodynamic model in patients with P. vivax malaria. K-^^ is 
the transit absorption rate constant (K^^ = (n + 1 )/mean absorption 
transit time), CL is the elimination clearance, is the apparent 
central volume of distribution, is the apparent peripheral volume 
of distribution, 0 is the intercompartment clearance, and Fis the 
relative oral bioavailability. 

nonlinear mixed-effects modeling. Substantial differences 
in matched venous and capillary blood concentrations have 
been reported in a previous clinical study.^^ A proportional 
venous-capillary transformation factor was, therefore, used 
to compensate for the difference between sampling matrices, 
resulting in 41% higher capillary concentrations compared 
with venous concentrations. 

The initial structural base models were parameterized as 
elimination clearance (CUF), intercompartment clearance(s) 
(Q/F), apparent volume of distribution of the central com- 
partment {VJF), apparent volume of distribution of the 
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Table 2 Population pharmacokinetic-pharmacodynamic parameters of the final model describing piperaquine in patients with Plasmodium vivax maiaria 



Venous estimates 



Capillary estimates 





Population 


95% CI" 


%CV for 


95% CI" 


Population 


Parameter estimates 


estimates' (%RSE ") 


for estimates 


IIV ("/oRSE") 


for IIV 


estimates* 


Pharmacokinetics parameters 












CUF(\/h) 


54.1 (5.56) 


48.0-59.7 


— 


— 


38.5 


VJF(\) 


2,770 (10.4) 


2,209-3,298 


— 


— 


1,970 


Q,/F (l/h) 


663(19.7) 


475-1 ,023 


138 (32.2) 


75.0-222 


471 


V„/F(l) 


10,600(10.7) 


8,549-12,962 


36.5 (40.2) 


18.8-51.7 


7,520 


Qj/F (l/h) 


72.7 (12.1) 


56.3-93.5 


— 


— 


51.7 


VJF(l) 


25,900(12.8) 


21 ,336-33,547 


45.2 (71.0) 


0.431-82.5 


18,500 


MTT (h) 


1.93 (9.16) 


1 .60-2.30 


54.5 (38.4) 


33.6-80.0 


1.93 


F(%) 


1 fixed 


— 


46.5 (18.9) 


36.0-55.1 


1 fixed 


TP 


1.41 (2.11) 


1 .35-1 .46 


— 


— 


— 




0.270 (5.59) 


0.239-0.299 










0.321 (5.79) 


0.283-0.354 








Pharmacodynamics parameters 












Bg^^ (relapses/year) 


8.94 (22.1) 


4.80-1 1 .7 








^s. (%) 


123 (23.7) 


61.8-165 








^sw (h) 


49.5 (54.3) 


35.7-80.4 








n 

^SHP 


4 fixed 










PC,„ (ng/ml) 


6.92 (15.5) 


6.15-11.3 








r 


9.28 (25.6) 


5.09-14.0 









CUF, elimination clearance; F, relative oral bioavailability; MTT, mean transit absorption time; VJF, apparent central volume of distribution; Q/F, intercompart- 
mental clearance(s); TP, venous-to-capillary transformation factor; VJF, apparent peripheral volume(s) of distribution; c,, residual error of piperaquine capillary 
concentrations; a,, residual error of piperaquine venous concentrations; 9^^^, baseline hazard; e^^, amplitude of surge function; 6^^,,, width of surge function; e^^p 
shape factor of surge function; PCj„, 50% protective venous concentration for 50% reduction in baseline hazard; r, shape factor of E^^^ function. 
■Computed population mean parameter estimates from NONMEM. Coefficient of variation (%CV) for interindividual variability is calculated as 100 x (EXP(mean 
estimate)-1)"^. 'Relative standard errors (RSE) are calculated as 100 x standard deviation/mean value from the nonparametric bootstrap method of the final 
pharmacokinetic model {n = 1 ,000) and pharmacodynamic model (n = 500). The 95% confidence intervals (95% CI) are displayed as the 2.5-97.5 percentiles of 
the bootstrap estimates. 



Table 3 Secondary parameters of the final piperaquine model in patients 



with Plasmodium vivax malaria 




Venous plasma 


Capillary plasma 


Cmax (ng/ml) 


162 (50.6-457) 


228 (71 .3-645) 


(h) 


3.0 (1 .6-9.3) 


3.0 (1.6-9.3) 


CL/F (l/h/kg) 


1.08 (1.01-1.66) 


0.770 (0.716-1.18) 


VJF{\Mg) 


802 (455-1,451) 


572 (324-1 ,033) 


',,2 (days) 


28.8 (15.5-58.3) 


28.9 (15.5-58.4) 


AUC^ (h X ng/ml) 


27,160 (9,443-78,470) 


38,250 (13,289-110,510) 


Day 7 (ng/ml) 


33.0 (11.4-95.2) 


46.5 (16.0-134) 



All values are given as median (range) unless otherwise specified. 
Secondary parameter estimates are calculated from the empirical Bayesian 
post hoc estimates. 

AUC^, area under the concentration-time curve from time 0 to infinity; 
C^^^, maximum concentration; CUF, elimination clearance; day 7, predicted 
concentration at day 7; T^^, time to maximum concentration; f^^, terminal 
elimination half-life; VJF, total volume of distribution. 

peripheral connpartment(s) (Vp/F), and absorption rate con- 
stant. Several absorption, distribution, and covariate mod- 
els were fitted to the data to construct the best performing 
model. Piperaquine pharmacokinetics in patients with P. vivax 
malaria were best described by a three-compartment distri- 
bution model (Figure 1). A three-compartment disposition 
model resulted in a significant improvement in model fit com- 
pared with a two-compartment model (difference in objec- 
tive function value (AOFV) of 108) with no additional benefit 
of an additional peripheral compartment (AOFV = 0.57). A 
transit-compartment model with a fixed number of transit 



compartments (n = 3, transit rate constant was set to be iden- 
tical to absorption rate constant) described the absorption 
phase better than all other absorption models (AOFV > 54.2). 
A relative bioavailability parameter (i.e., fixed to 100% for the 
population) was implemented to allow quantification of the 
interindividual variability in the absorption of piperaquine and 
resulted in a significant improvement in the model fit (AOFV = 
81.8). 

Body weight, incorporated as a fixed allometric function, 
improved the model fit significantly (AOFV = 31.0) and resulted 
in better precision and decreased interindividual variability 
in clearance and volume parameters. Body temperature on 
admission on VJF(AOF\/ = 1 1 .5) and age on Q/F (AOFV = 
5.46) were selected in the forward covariate selection but 
could not be retained in the backward deletion. The estimated 
interindividual variability was small for the transformation fac- 
tor, CUF, QJF, and VJF, and could, therefore, be removed 
with no significant impact in the final model (AOFV = 2.40). 
Final pharmacokinetic parameter estimates and secondary 
parameter estimates are reported in Tables 2 and 3. 

The mean parameter estimates from the simultane- 
ous model with both venous and capillary data were also 
comparable with estimates when capillary and venous 
data were modeled separately. The final pharmacokinetic 
model displayed satisfactory goodness-of-fit diagnostics 
(Supplementary Figure SI online) with a small epsilon- 
shrinkage of 15.0%. However, a relatively high ri-shrinkage 
(up to 43.5%) could be seen for certain parameters in the 
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Figure 2 Prediction-corrected visual predictive cliecl< of the final 
model describing piperaquine pharmacokinetics in P. wVax malaria. 
Open circles represent the observed piperaquine concentrations. 
Solid black line represents the 50th percentile of the observations, 
and dash lines represent the 5th and 95th percentiles of the 
observations. Gray areas represent the 95% confidence intervals of 
the simulated 5th, 50th, and 95th percentiles from 1 ,000 simulations. 
The inset shows a prediction-corrected visual predictive check for 
the first 3 days of piperaquine treatment. 

final model because of the sparse data. The numerical and 
prediction-corrected visual predictive checks resulted in 
4.6% (95% CI: 2.7-8.1%) and 3.7% (95% CI: 2.5-8.2%) of 
piperaquine observations below and above the simulated 
90% prediction interval, respectively (Figure 2). 

Pharmacodynamics of piperaquine 

Piperaquine inhibits asexual parasite multiplication but has 
no effect on liver-stage parasites. The influence of pipera- 
quine exposure on the risk of recurrent P. vivax malaria was 
investigated in 62 patients for whom both pharmacokinetic 
and pharmacodynamic data were available. However, effi- 
cacy data were available from all the 241 patients included 
in the study. Pharmacodynamic parameters were evaluated 
using an interval-censoring time-to-event model, imple- 
mented with the Laplacian estimation method. The phar- 
macokinetic parameter and variability estimates were fixed 
to that of the final pharmacokinetic model and used in the 
pharmacokinetic-pharmacodynamic model. P. vivax malaria 
in tropical regions commonly displays frequent relapses 
which emerge at three-week intervals when the infections 
are treated with rapidly eliminated antimalarials.^^ A multiple 
surge function was, therefore, implemented to characterize 
this periodicity in the risk of recurrent malaria. The constant 
hazard function with a periodically increased risk of relapses 
described the data well. Piperaquine had a significant (AGFV 
= 20.3) protective effect on the risk of recurrent malaria infec- 
tions with a required in vivo venous piperaquine concentra- 
tion of 6.92ng/ml for a 50% relative reduction in baseline 
hazard (PC ). Sex was selected as a significant covariate 




Time (days) 




30 40 
Time (days) 

Figure 3 Visual predictive check of the final time-to-event model 
describing piperaquine pharmacodynamics in P. wVax malaria. Gray 
area represents the 95% prediction interval of the time to recurrent 
infections from 500 simulations. Solid line represents the observed 
Kaplan-Meier plot for (a) patients in the pharmacokinetic and 
pharmacodynamic cohort (62 patients. Table 1 ) and (b) all patients 
in the efficacy study (241 patients. Table 1). 

(P < 0.05) on baseline hazard in the forward selection step, 
but it was removed during the backward elimination resulting 
in a covariate-free pharmacodynamic model. Final pharma- 
codynamic parameters were well estimated with high preci- 
sion (Table 2) and simulation-based diagnostics showed 
good agreement between simulated and observed recurrent 
malaria infections (Figure 3). 

The individual patient parasite burden (Table 1) corre- 
sponded approximately to a total of 10" parasites in a typi- 
cal adult at admission. This assumes that P. vivax does not 
sequester substantially. DHA-PQ is administered for 3 con- 
secutive days, covering two asexual cycles. Assuming a 
parasiticidal effect (parasite reduction ratio) of a 10,000-fold 
parasite reduction per cycle by dihydroartemisinin (without 
any additional effect by piperaquine) results in a 100,000,000- 
fold reduction in total parasite biomass during the first two 
cycles.^" At an initial parasite burden of 10" parasites, this 
would leave -1,000 parasites to be eliminated by residual 
piperaquine concentrations to prevent recrudescent malaria. 
Piperaquine alone has not been adequately assessed in 
P. vivax malaria, but if it is similar to chloroquine, it should 
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Figure 4 Relationship between piperaquine treatment and risk 
of relapsing P. vivax malaria infections, (a) Mean hazard rate of 
relapsing P. vivax malaria versus time in the presence (solid line) 
and absence (dashed line) of a standard piperaquine treatment 
(primary y-axis). The frequency of observed recurrent P. vivax 
infections for all patients in the efficacy study (241 patients, Table 
1) is displayed on the secondary y-axis. The hazard of the first 
relapse (3 weeks) is completely suppressed and the hazard of 
the second relapse (6 weeks) is almost completely suppressed 
after piperaquine treatment. The model-predicted hazard rate 
after piperaquine treatment is in agreement with the observed 
frequency of observed recurrent P. wVax infections, (b) Simulated 
day-7 venous piperaquine plasma concentrations and chance of 
remaining malaria free for 30 days. Open gray circles represent 
simulated patients. Solid black line represents the 50th percentile 
and dash lines represent the 5th and 95th percentiles of the 
prediction interval. The dashed vertical line indicates the predicted 
day-7 concentration that results in 99% chance of remaining malaria 
free for 30 days for a typical patient. 

provide parasite reduction ratios of -100 to 1,000 per cycle, 
so most patients should be parasite free 6-8 days post-treat- 
ment if piperaquine levels are sufficient to sustain a maximum 
parasiticidal effect for this time. If the parasite burden is not 
eliminated completely, then these remaining parasites can 
multiply and cause a recrudescent infection. However, a high 
cure rate of the asexual infection is likely in this study as con- 
centrations of piperaquine after the second asexual cycle (i.e., 
> 4 days) are highly likely to be able to eliminate ~1 ,000 para- 
sites. In the first relapse of tropical P. wVax malaria (unaffected 
by drugs and assuming 1-10 schizonts and a subsequent 
parasite multiplication rate of -10 per cycle), the hypnozoite- 
derived hepatic schizonts must liberate 13-17 days after the 
onset of the primary illness to produce a symptomatic infec- 
tion at day 21 . However, if residual piperaquine concentrations 



are sufficiently high, they could suppress and eliminate the 
asexual parasites derived from 10,000 to 100,000 released 
merozoites and thereby prevent the first relapse completely. 
If residual piperaquine concentrations are not sufficiently high 
for complete elimination, they could cause an initial suppres- 
sion, resulting in a delayed symptomatic first relapse. Only 
14% (15 out of 109) of recurrent P. vivax infections occurred 
before day 39, which suggests that DHA-PQ treatment pre- 
vents a large proportion of the first relapses (Figure 4). The 
majority of infections (64%, 70 out of 109) occurred between 
day 42 and 56, suggesting that these are mainly delayed 
second relapses (Figure 4). There will be some contribution 
from reinfections, but as entomological inoculation rates are 
probably well below 0.5/person/year for P. vivax in this area, 
this contribution is likely to be small. 

DISCUSSION 

The fixed combination of DHA-PQ is a highly efficacious 
antimalarial drug with an excellent safety profile. Antimalar- 
ial treatment with DHA-PQ resulted in lower cumulative risk 
(54.9%) of recurrent P. vivax malaria compared with chloro- 
quine (79.1%) in Thailand.^'' This difference could be a result 
of the very potent artemisinin derivative which generally has 
a rapid onset of action against P. wVax malaria, but it is likely 
to result mainly from increasing chloroquine resistance and 
hence the superiority of piperaquine. In fully sensitive P. vivax 
malaria, recrudescence rates are close to zero with chloro- 
quine treatment. Previous studies in this region have shown 
high relapse rates of P. vivax within one month of treatment 
with a short acting artemisinin derivative. Thus, both slowly 
eliminated drugs prevent the first relapse. Whether the lower 
cumulative risk of recurrent P. vivax malaria, following DHA- 
PQ, results from a higher cure rate of the primary infection or 
a more effective suppression of the first or the second relapse 
has been unclear. As relapses can arise from parasites which 
are either genetically similar or genetically distinct to those 
causing the primary infection, it is not possible to distinguish 
with certainty between relapse, recrudescence, or reinfection 
in patients with P. vivax malaria who remain in the endemic 
area. Recrudescence can occur many weeks after treatment 
with slowly eliminated drugs, but the pharmacokinetic-phar- 
macodynamic data modeled in this study suggest that only a 
negligible proportion (<2%) of the recurrent P. wVax infections 
could be recrudescent. 

Pharmacokinetics of piperaquine 

The population pharmacokinetic properties of piperaquine 
have been described previously in patients with uncom- 
plicated P. falciparum malaria. This is the first study 
reporting the population pharmacokinetics of piperaquine in 
patients with P. vivax malaria. 

The structural three-compartment model described in this 
study was similar to that reported recently in young children 
and in pregnant and nonpregnant women with P. falciparum 
malaria. ^^'^^'2'' Simplified two-compartment structural models 
in children and adult patients have also been reported reflect- 
ing differences in follow-period and/or sampling strategies. 
A transit-compartment absorption model allowed more 
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flexibility in the absorption phase and provided a clear advan- 
tage over other absorption models, although there were not 
enough data in the absorption phase to estimate both K^^ 
and separately. Large interindividual variability was seen 
in the absorption of piperaquine as in previous studiesJ^'^^'^^ 
Data collected here were too sparse to allow an estimation 
of the variability between dose occasions, which is likely to 
have inflated the inter-individual variability. Indeed interin- 
dividual variability was high in this study with values up to 
138%. Venous and capillary plasma concentrations were 
successfully modeled simultaneously with a transforma- 
tion factor of 1.41 (relative standard error of 2.11%). This 
factor was similar to that from a linear regression between 
observed venous and capillary plasma concentrations (slope 
= 1.29), which further supports the appropriateness of the 
simultaneous model. However, the exact mechanism of this 
matrix-dependent difference cannot be elucidated from the 
data collected in this study. The advantage of a simultaneous 
modeling approach is that the model can be used to simulate 
drug exposures from any sampling technique and therefore 
enables literature comparisons. 

Body weight was the only significant covariate in this 
model. This has also been reported in previous studies^^-^'* 
and was not unexpected, considering the strong biological 
prior of body weight as a covariate on pharmacokinetic dis- 
position and elimination parameters. All pharmacokinetic 
parameters were estimated with good precision (relative 
standard error < 20%) and the simulation-based diagnostics 
indicated excellent predictive performance of the final phar- 
macokinetic model. 

Piperaquine CUF in this P. vivax study was 1 .08 l/h/kg, 
total volume of distribution {VJF} was 802 I/kg, and f,,^ 
28.8 days, which were comparable with previously pub- 
lished reports in patients with an uncomplicated P. falciparum 
malaria (mean (range) CUF oi 1.18 (0.90-1.4) l/h/kg, VJF 
of 789 (574-877) I/kg, and f,,^ of 23.7 (18-28) days).^^^-^''-""'^'* 

Pharmacodynamics of piperaquine 

The observed high relapse rate of P. vivax malaria was 
expected as DHA-PQ eliminates only blood-stage parasites 
and does not provide cure of the latent liver-stage parasites. 
It was expected that piperaquine would delay the time to 
relapse, as does chloroquine, and this was supported by 
the pharmacodynamic results. Tropical frequent relapsing 
strains of P. vivax typically relapse at 3-4-week intervals, 
and the median time to relapse approximates three weeks 
(21 days) in Thailand after treatment with rapidly eliminated 
antimalarial drugs (such as artesunate, quinine, and halo- 
fantrine).^" In an area of much higher transmission in Papua, 
Indonesia, there were significantly more early recurrences 
of P. vivax malaria following artemether-lumefantrine treat- 
ment of P. vivax malaria compared with DHA-PQ (cumula- 
tive risk of 38 vs. 10% at day 42)^^ reflecting the failure of 
the more rapidly eliminated lumefantrine to suppress the 
first relapse. However, with longer follow-up in a high trans- 
mission area, the cumulative reinfection rates converge as 
everyone is reinfected eventually. A critical and unanswered 
question has been whether the slowly eliminated drugs pre- 
vented, or simply delayed, the relapse. Explicitly, prevention 
would mean that the relapse emerged but was eliminated 



by the residual drug levels, and thus, the first observed 
recurrence is in fact the second relapse. Delay would mean 
that the residual drug levels reduced multiplication of the 
relapse, so it reached patency after a substantial delay. In 
the former explanation, overall relapse numbers might be 
reduced as this would exhaust the liver burden of hypnozo- 
ites more rapidly. A pharmacometric time-to-event approach 
explaining relapsing P. vivax malaria, first presented for 
amodiaquine treatment in pregnant women, described 
the data observed in this study well. An interval-censoring 
time-to-event model was successfully used since the exact 
times of relapses were unknown. The risk of relapsing P. 
vivax malaria in patients was explained by a constant base- 
line hazard (i.e., 8.94 relapses/year) with a 123% increased 
risk every third week. The pharmacodynamic model was 
significantly improved by adding a protective drug effect, 
indicating that piperaquine cured the primary blood stage 
P. vivax infection and fully suppressed a large proportion 
of the first relapses that would have occurred three weeks 
after the initial infection. The pharmacodynamic model also 
indicated that the observed high recurrence rate is likely to 
be explained mainly by delayed second relapses, although 
it cannot be excluded that some recurrences are reinfec- 
tions. This can only be mechanistically evaluated by an 
individual parasite-data-driven pharmacokinetic-pharma- 
codynamic model. Simulations using the final pharmaco- 
kinetic-pharmacodynamic model (Figure 4) demonstrated 
that minimum day-7 piperaquine venous and capillary con- 
centrations of 27 and 38ng/ml, respectively, were needed to 
suppress the risk of relapse for 30 days (i.e., for 99% of sim- 
ulated patients). This is in close agreement with a previously 
defined day-7 threshold venous concentration of 30ng/ml in 
patients with P. falciparum malaria to minimize the risk of 
recrudescent infections.^" 

In conclusion, the population pharmacokinetic proper- 
ties of piperaquine in patients with P. vivax malaria were 
described successfully by modeling venous and capillary 
plasma concentrations simultaneously. Body weight was 
the only significant covariate and resulted in increasing 
clearance and volumes with increasing body weight. Times 
to recurrent P. vivax infections were successfully modeled 
with a time-to-event approach and resulted in a significantly 
delayed time to recurrent infections during the post-treat- 
ment phase of piperaquine treatment. This was explained 
by prevention of the first P. vivax malaria relapse by residual 
piperaquine concentrations. Piperaquine is a good candi- 
date for the treatment of P. vivax malaria and the model- 
ing conducted here demonstrated the added benefit of 
reduced morbidity compared with other more rapidly elimi- 
nated antimalarial drugs. 

IVIETHODS 
Ethics approval. 

Ethics approval was granted by the Faculty of Tropical 
Medicine Ethics Committee, Mahidol University, Bangkok, 
Thailand, and the Oxford Tropical Research Ethics Commit- 
tee, Oxford University, UK. This study was registered in the 
ISRCTN Register (ISRCTN87827353). Details on the clini- 
cal study design and outcomes have been presented in full 
elsewhere. 
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Drug regimen and blood collection. 

All patients received a standard regimen of DHA-PQ. Patients 
who vomited the dose were excluded from the pharmacoi<i- 
netic study. A full or half replacement dose was administered 
if vomiting occurred within 30min of or between 30min to 1 h 
after administration, respectively. 

Venous blood samples (3 ml) and capillary blood samples 
(200 ^Jl) were drawn from each patient for piperaquine plasma 
measurements. Samples were collected randomly over 69 
days (4-7 venous samples and 1-11 capillary samples) 
from 62 patients and an additional venous sample (5 ml) was 
drawn from these patients at the time of recurrent P. vivax 
malaria. This random sampling allowed of an adequate cov- 
erage of the entire concentration-time profile. An additional 
54 patients provided a single plasma sample only at the time 
of recurrent malaria. Patients were excluded from the study 
after recurrent malaria. Plasma samples were shipped on dry 
ice to the Department of Clinical Pharmacology, IVIahidol- 
Oxford Tropical Medicine Research Unit, Thailand, for drug 
quantification.^^ 

Population modeling. 

Piperaquine pharmacol<inetics and pharmacodynamics 
were evaluated using nonlinear mixed-effects modeling 
(NONIVIEM version VII; Icon Development Solutions, Ellicott 
City, MD). Postprocessing and automation were performed 
using Pearl-Speal<s-NONMEI\/l version 3.5.3, Census ver- 
sion 1.2b2,^^ and Xpose version 4.0^^ in the programming 
language R version 2.13.2 (The R Foundation for Statistical 
Computing). Measurements below the lower limit of quanti- 
fication were less than 0.5% of the total data and therefore 
omitted. OFV, calculated by NONMEM as minus twice the 
log-likelihood up to constant, was used for model selection 
during the model-building process. A difference in the OFV 
(AOFV) of >3.84 was considered significant (P< 0.05) when 
comparing two nested models with 1 degree of freedom 
difference. Full details of the pharmacol<inetic and pharma- 
codynamic modeling methodology can be found in supple- 
mentary material. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

^ Plasmodium vivax malaria causes substantial 
morbidity in transmission areas as a result of 
multiple relapses. P. vivax resistance to chloro- 
quine is increasing in many parts of \he world. 
Diliydroartemisnin-piperaquine is a promising 
candidate for the treatment of P. vivax malaria. 



WHAT QUESTION DID THIS STUDY ADDRESS? 

^ Tine aim of this study was to evaluate the popula- 
tion pharmacokinetics and pharmacodynamics 
of piperaquine in patients with P. vivax malaria 
in Tliailand and to determine whether residual 
piperaquine prevents or delays the emergence 
of P. vivax relapse. 

WHATTHIS STUDY ADDS TO OUR KNOWLEDGE 

^ Piperaquine pharmacokinetics and the biology 
of P. vivax malaria were successfully modeled. 
A large proportion of the first relapses were 
suppressed completely by residual piperaquine 
concentrations, and recurrences resulted main- 
ly from the second or the third relapse or from 
reinfection. 



HOW THIS MIGHT CHANGE CLINICAL 
PHARIVIACOLOGY AND THERAPEUTICS 

^ Dihydroartemisinin-piperaquine is a promising 
candidate for the treatment of P. vivax malaria 
and the modeling conducted here demonstrat- 
ed the added benefit of reduced morbidity com- 
pared with other short-lived antimalarial drugs. 
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